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[57] ABSTRACT 

The size of a digitally-stored image is modified as de- 
sired by expanding the source image data and adjusting 
the input to the output raster display device accord- 
ingly. 

12 Claims, 3 Drawing Sheets 
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4,080,634 or 4,067,021 or a CRT exposure tube. Simi- 

METHOD AND APPARATUS FOR PRODUCING A larly other printers employing a thermal head of the 

SPECIFIED FORMAT OUTPUT IMAGE FROM AN type disclosed in U.S. Pat. No. 3,951,247 may be used 

ARBITRARY FORMAT SOURCE IMAGE and the like. , 

FIFT n OF THF TKVPWTTnW 5 A ***** P r0duCCd ^ SUch device ' is K ener " 

hltLU Uh rHE INVENTION atcd by writing with a moving spot of controlled (mod- 

This invention relates to a method and associated ulated) energy on a sensitive film sheet, point by point, 

apparatus for altering the output displayed size of digi- along sequential raster lines, spots of variable density 

tally stored image and more particularly for electronic which correspond to the original pixel information that 

cally altering the size of a stored digital image so as to 10 was extracted from an original image and stored in a 

reproduce it using a raster display in a predetermined storage medium. 

size and format This type of output display device is well known in 

BACKGROUND OF THE ART ? C ne ? d n ? 1 * cIaborated u P° n hcre > but will 

be referred to hereinafter as a recorder, (regardless of 

Digital image storage is a process whereby an image 15 whether the raster is generated using electrical, me- 

is stored in a storage medium as a plurality of bits repre- chanical or combination methods, and without regard 

senting pixels (picture elements) of a given intensity. to the nature of the writing source, provided such 

Digital storage offers a number of distinct advantages to source is amenable to some form of output modulation 

image storage and transmission, image processing, control which alters the resulting density of the writing 

image retrieval etc. over the more traditional hard 20 0 n an output film sheet) ■ 

copy. From the point of ready access, a computerized To reproduce an original picture in exactly the same 

nl^fwn'^r f M ^ ^ ^ n *T !° ^ * ™<**" * SDOt size the 

ulatc than paper folders containing film sheets. Image size as that of the scanner first used to acquire the 

ttansmiss.on over telephone lines to remote locations stored image puels. Thus a 3 inch by 3 inch rectangular 

S^r^t 11 15 n ° l ^™ ""P™* r th f 25 image acquired with a scanner having a resolution^ 

digital image capture, conversion and storage, finds im «i*aw*r ™v i Z 

great utility in applications where there is need for such JJ ^^ < ^^?? Wf? 

advantages g 0 g X fa 01120111 * 11 ) *nd Y (vertical) directions 

Applications are found in fields otherwise as diverse ^** tcnsive ™? h ^"ft b ° rdcrs ' reproduced 
as medical diagnostic radiology and the publishing in- 30 cxac ^ ,n r a ? \$ 3 «ch space using a recorder 
dustry, where electronic page composition is becoming capable of writing 300 spots and 300 lines per inch. The 
increasingly common. Often in such applications there fittm * L°J the ° ngmal pi . CtUrc m 016 dlotted 8pace wi!1 be 
is need to alter the original size of a given image data to straightforward, one pixel to one spot and no data will 
fit within prescribed space limitations. These limitations b v e ]os } or created. However, if, for example, 
are imposed by a number of extraneous factors which 35 the onginal image is composed of 900 pixels X 900 raster 
are usually the result of physical constraints in the out- lmes » which must f,tted a 2x2 inch space using a 
put product, e.g., a sheet of film or paper on which one recorder having a resolution of 1000 spots/inch and 
or more images are displayed. For instance, in a medical 1000 raster lines/inch, the resulting output image size in 
diagnostic application, it may be advantageous to dis- tcrms of s P ots 211(3 lines equals 2000 x 2000 spots for 
play a number of different radiograms, side by side on a 40 which we to provide intensity values. Therefore there is 
single radiological film for a physician to examine. Simi- insufficient data from the original to fill the excess spots, 
larly, a number of distinct images and corresponding ° ne solution used in the prior art is to generate ficti- 
text may have to be fitted on a sheet of paper forming a. tious pixel values from those on hand, The simplest way 
magazine page. In both instances, it is rare that the m tn * s instance is to double the original pixels by repli- 
original images all happened to be the exact size that 45 eating each at least once so that we now have an origi- 
will fit in the available space. In most cases the stored nal image composed of 1800x1800 pixels and lines, 
image size will have to be altered to a new size commen- This, however, only works for magnifications that are 
surate with the space it is to occupy in the display me- exact multiples of the original. Another approach is to 
dium. try to predict values to generate imaginary pixels to fill 
One possible solution, albeit cumbersome and slow, is 50 the voids as shown in U.S. Pat. No. 4,595,958. How- 
to generate an image that is an exact duplicate of the ever, there is a problem in the creation and allocation of 
original using the digital data, and then project this pixels for the partial magnification, i.e., the need to 
image using variable magnification lens systems onto a magnify the image 2000/900=2.222 times. The genera- 
portion of a photosensitive film sheet. By varying the tion of data for, and the distribution of, the extra 
magnification the projected image can be made to fill 55 200x200 spots in the output image end up creating 
the allocated space. The process is then repeated for artifacts and/or image fuzziness in the final product, 
each of the images and sizes used until the film sheet is U.S. Pat. No. 4,080,634 briefly discusses a method to 
Tilled. scale original analog images for digital recording. Ac- 
One would rather, however, employ electronic cording to its teachings, the scanning (reading) rate of 
means which operate directly and rapidly on the digital 60 an original image may be changed so that by the time 
data to provide the displayed image data in proper for- the spot reaches the end of each line, a preselected 
mat so that the displayed image fits exactly in the alio- number of samples have been taken for that line. At the 
cated space. That is one would rather have direct, elec- same time, the rate of advancement of the image under 
tronic image magnification for image fitting and the the scanning line is modified so that the vertical scan 
ability to do so with or without altering the original 65 distance compensates for the horizontal reading rate to 
image aspect ratio as desired. The image data could then maintain the original image aspect ratio, 
be used to drive an output device which may be a laser The problem with such a process as applied to re-, 
printer scanner such as described in U.S. Pat. Nos. cording, is that the residence time of the spot per unit 
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area changes, which in turn changes the exposure of the 
writing or recording medium. For small changes in 
magnification, this is not a problem, but for changes 
such as 2, 4, 8, or 10 times magnification, one starts 
encountering exposure problems, due to reciprocity 
failure and incorrect exposure in photosensitive sys- 
tems, or thermal delay in thermally activated systems. 
The lines in the vertical direction become separated and 
the spaces between lines start becoming visible. There is 
still need for an effective image altering system capable 
of fitting an arbitrary digitally stored image retrieved 
from a storage medium into an arbitrary fixed allocated 
area in a display medium, which will not result in loss of 
information or introduction of artifacts or distortion. 

SUMMARY OF THE INVENTION 

The aforementioned difficulties may be resolved 
through a process for use in a raster scan type image 
reproduction system for fitting a source image com- 
prised of source image data representing an arbitrary 20 
number of pixels and an arbitrary aspect ratio in a speci- 
fied output image having a desired size and aspect ratio 
comprising the steps of: 

(a) determining a desired output image size in at least 
one dimension, 25 

(b) determining a source image size in the same at 
least one dimension, and comparing said source 
image size to said output image size to obtain a ratio 
R in that at least one dimension 

(c) calculating a first integer I adjacent, preferably 30 
next higher, to the ratio R 

(d) digitally expanding the source image by expand- 
ing the source image data by a factor equal to the 
calculated integer I to generate an expanded output 
image data 

(e) calculating a sampling factor f where f=I/R, and 
(0 Sampling the expanded output data at a sampling 

rate controlled by the sampling factor f and use the 
sampled output data to drive an output raster dis- 
play device, such as a recorder. 
Still within the scope of the present invention, the 
raster scan system of the above process may scan a 
surface at a specified scan rate to generate an X-Y ras- 
ter, and the factor f=I/R may be used used to control 



source image data represented by an arbitrary number 
of pixels and an arbitrary aspect ratio, the data residing 
in a storage means, and for generating data for repro- 
ducing in said recorder said image in a specified output 
image format having a desired size and aspect ratio 
comprising: 

(a) a central data processing means including means 
for accessing a preselected image data in a storage 
medium; 

(b) means also in said central processing means for 
calculating: 

i. a ratio R where R= output image size/source 
image size; 

ii. a first integer I adjacent, preferably next higher, 
to the ratio R; and 

iii. a quotient of the ratio R to the calculated inte- 
ger I to obtain a sampling factor f where f=I/R; 

(c) means controlled by said central processing 
means, for retrieving said preselected image data 
and for delivering said data to a data interpolating 
and storing means adapted to digitally expand the 
source image data by a factor equal to the calcu- 
lated integer I to generate expanded output image 
data and to store at least a portion of such data; and 

(d) means associated with the central data processing 
means to apply the sampling factor f to a sampling 
rate control means to adjust the sampling rate of 
the expanded output data generated by the interpo- 
lator, and to supply said data to the recorder, and 
to also apply the sampling factor f to the control 
means for adjusting the vertical scanning rate of 
the recorder. 



35 



such raster scan rate in both X and Y directions. Where 45 and 



BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will next be described with reference 
to the attached drawings in which: 

FIG. 1 shows a greatly enlarged portion of a raster 
image produced by a raster scan record, 
FIG. 2 shows a greatly enlarged portion of an image 
40 as it is scanned by a readout scanner. 

FIG. 3 shows a table used to store pixel image infor- 
mation. 

FIG. 4 represents in schematic format an apparatus 
for implementing the process of the present invention, 



the output data is used to modulate a writing beam 
intensity, the intensity of the beam may also be adjusted 
to compensate for changes in exposure resulting from a 
corresponding change in the sampling rate. 

The process in accordance with the present invention 50 
may further comprise implementing the digital expand- 
ing step using any one of the following methods: 

(a) Replication, or 

(b) By-linear interpolation, or 

(c) Cubic spline interpolation, or 55 
(e) Two dimensional kernel convolution interpola- 
tion. 

The process of this invention may be advantageously 
employed in instances where the desired output image 



FIG. 5 shows a film sheet on which 9 (nine) images 
are displayed, arranged in accordance with the teach- 
ings of this invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The invention will next be described with reference 
to the drawings which are used herein for illustration 
and as such are not intended to limit the invention to the 
specific embodiments depicted. 

Referring first to FIG. 1 there is shown in schematic 
representation a greatly enlarged portion of an image 
reproduced using a raster scan recorder. The image is 
reproduced using a focussed beam of energy traveling 



size is determined based on optimum distribution of a 60 across an imaging medium which upon exposure to the 



given number of source, images on a preselected size 
film sheet and more particularly, wherein the source 
images represent medical diagnostic images. 

The present invention further comprises, associated 
apparatus for use with a raster scan type image repro- 
duction recorder having control means to adjust a scan- 
ning rate in at least one dimension, the equipment 
adapted for accessing a source image comprised of 



beam energy produces a visible image. The beam inten- 
sity is modulated to expose the imaging medium with 
variable energy levels. The image created on the me- 
dium may be viewed as comprising a series of spots 64 
65 laid along a raster line 62 extending along one side 67 of 
a generally rectangular sheet 65. The spots are written 
on the sheet along the line 62 from left to right in a 
direction defined arbitrarily as the X or horizontal scan 
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direction. The full surface of the sheet 65 may be FIG. 3 shows how the output of the scanner/digitizer 
viewed as covered with spots written along raster lines may be stored in an image storage system. For each 
that arc placed adjacent each other in a direction 68, pixel there is an address, e.g., ^5", and an optical den- 
generally perpendicular to the X direction arbitrarily sity value. This optical density value is represented by a 
defined as the Y or vertical scan direction. Individual 5 digital number which in an 8 bit binary system ranges 
spots are of course only detectable in a limiting case in from 0 to 255. The address and the density value may be 
which the modulation frequency and pattern of the stored in a table as illustrated, for later access. The 
beam intensity and the scanning rate have been selected numbers of course could also be simply stored in plain 
to generate image spots. More typically in actual imag- serial form, with pointers to indicate end of lines, etc. 
ing situations the spots tend to blend together as the 10 The storage medium may be either a long term (such as 
beam scans the imaging or writing medium. However tape, magnetic disk or optical disk) or volatile or perma- 
the individual spot concept is maintained herein as an nent electronic memory. 

aid in explaining the problem and solution offered by FIG. 4 shows an arrangement of components includ- 
the present invention. For the same purpose, the fact ing both hardware and software for practicing the pres- 
that the spot intensity is not constant across a spot diam- 15 ent invention. The system may comprise a Central Pro- 
cter but falls of from the center outwardly usually in a cessing Unit ("CPU 12") which may be a personal com- 
Gaussian distribution manner, is also ignored, and the puter programmed to perform the necessary calcula- 
typical practice of taking an arbitral^ point in the inten- tions and provide the needed output signals as will be 
sity curve, such as J of max., is used, creating an admit- explained below. 

tedly artificial but conceptually useful spot with distinct 20 At least two information items are inputed, the output 

limits. All these concepts are well known in the art and system parameter limitations, i.e., the system resolution 

are used to facilitate visualization of the claimed pro- in terms of spots per inch, and the output overall size in 

cess, rather than limiting such process. terms of total spots in each orthogonal (X-Y) direction. 

By studying the spot pattern shown in FIG. 1 it is Line 14 indicates this input. The second input comprises 

seen that a recorder's resolution, i.e., the smallest indi- 25 the image size to be reproduced again in terms of total 

yidually distinguishable area of information is substan- pixels in the X and Y directions. This input is repre- 

tially the same as the smallest writing spot possible with sented by an line 16 in FIG. 4. ■ 

that recorder, assuming that the modulator bandwidth The source of information regarding the output size 

is high enough so that a single spot written along the of the reproduced image is not shown in FIG. 4 since it 

scan line will appear round on the exposed medium. It 30 is not essential for the practice of this invention. In a 

has thus become customary to identify a recorder's preferred application, where the reproduced images 

performance in terms of number of raster lines per inch, represent radiological images and where it is desirable 

or if the spot is maintained symmetrical, in spots per to display a number of such images on a single film 

inch writing capability, rather than in terms of minimal sheet, the following procedure is followed to determine 

spot diameter size. For instance, a recorder may be 35 the output image size. 

given as having a resolution of 1000 (one thousand) Refer next to FIG. 5 in which an output film sheet is 

spots per inch and being capable of handling sheets up shown in reduced size. The film sheet is first identified 

to 10 (ten) inches wide. In describing this invention we in terms of physical size, e.g., 11 by 13.5 inches. The 

will use spots and spots per inch to identify physical resolution of the recorder is noted, e.g., 300 lines per 

image sizes. 40 inch (and 300 spots per inch). Any display limitations 

Similarly to image reproducing recorders, the size of are next considered. For instance, in order to permit 

a digitally captured image may be conveniently repre- clamping the film on a light source, a margin of three 

sented in terms of the number of picture elements avail- quarter inches (0.75") is allowed at the top. This corre- 

able for a given image size. Thus in the case where an sponds to 225 raster lines (75 inches ♦ 300 lines per inch) 

original image 55 shown in FIG. 2 is scanned with a 45 which are measured off to be left blank at the top of the 

moving spot 50 along a plurality of raster lines 52 and sheet, 

the output digitized, the result is a series of numbers The balance 3300 horizontal spots by 3825 vertical 

corresponding to the smallest picture element distin- raster lines define the available image area for the dis- 

guished by the scanning spot. This picture element or play of images. In this illustration, it is desired to display 

pixel 56 is at the scanning system limit similar to the 50 nine radiograms each with an aspect ratio of 0.84375. 

scanning spot size of a recorder, and as in the case of The original aspect ratio is to be maintained, i.e., the 

recorders, reading or analyze scanners/digitizers are magnification is to be the same in both X and Y dimen- 

also identified by their resolution capabilities in terms of sions. Some margins and image separation are desirable, 

pixels per inch and total reading width capabilities. Using the technique to be described in detail below, it is 

In practicing this invention, the source of the date 55 determined that the images will fill the sheet in an opti- 

representing an image need not have been generated mal way that is achieving maximum utilization of the 

with a scanner/digitizer and a preexisting image. The available film sheet surface if each image is displayed as 

same concepts apply in cases where the image is gener- 3.375 inches horizontally by 4.0 inches vertically. This 

ated in raster format directly on a detector and immedi- corresponds to each image occupying the space of 

ately converted into digital format, as may be the case in 60 1012.5 spots horizontally and 1200 raster lines verti- 

certain forms of medical radiology, such as tomography cally. 

and the like. Also, the image data may have been gener- The above information is inputed through input 16 to 

ated completely in a computer using for instance a the CPU 12. Of course if this calculation is done in the 

drafting program. The same principles and definitions same computer as is likely the case, there is no need for 

apply, each image being defined as being a particular 65 an external input of the result which can preferably be 

size the size given in terms of number of pixels in two internally supplied to the CPU for further processing, 

generally orthogonal dimensions. The source of the Other information is also provided to the CPU through 

image data is not important. additional inputs not shown here, such as data identify- 
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mg the image to be retrieved for display, the position of nine adjacent address values all with the same magni- 

the displayed image relative to the other images on the tude, three in adjacent horizontal addresses and the 

film sheet, patient identification text data, or other perti- same three in each of three sequential lines in effect 

ne J" date to be also displayed, etc multiplying the incoming data by an integral factor in 

Once this ^ information is ^ provided o the CPU, the 5 both the horizontal and vertical dimensions, 

image is retrieved from the image data m some memory As a resuJt , the temporary buffer now hoWs „ ^ 

interpolation step occurs, after which data is temporar- * e 'T^/P 0 ? StlJ1 " 0t fit thc desired ° tttput 

Uy stored in a buffer for use in displaying the image. The *™ 0f 1012,5 honzontal s P°ts °y 1200 vertical raster 

buffer may hold a completely expanded image, or in a 5 s * 

preferred embodiment in order to limit memory needs, In a ty P lcal recorder arrangement, a radiation beam 
it may contain a limited number of lines at a time, and 15 4 * "^tted from * radiation source 38 is directed 

data processing, including input and output may occur through a modulator 40 as a modulated beam 44 to 

simultaneously. The interpolator may further include impinge on a rotating mirror 46 driven by driving 

buffer storage means to receive and hold data from means 48 and therefrom onto the surface of moving film 

memory 20 prior to the interpolation step, in addition to sneet ^ which it scans along a line generally transverse 
the buffer means for holding data after interpolation, jo to the direction of movement of the film. The film sheet 

CPU 12 performs a number of operations in addition nas ; a surface sensitive to the radiation of the beam, 

to selecting the images to be displayed. First the size of Driving means for moving the film sheet in a path past 

the original image selected is determined based on the the scanning radiation beam 44 typically comprise rol- 

stored image pixel data. In the present illustration, we lers 34 and a driving motor 32. 
will assume that the original data was captured in a 25 In addition to calculating integer I CPU 12 is further 

typical TV monitor using a frame grabber, and that the programmed or otherwise adapted to calculate a second 

image comprises 432 horizontal pixels by 512 vertical factor f=I/R. This factor is used to control the rate at 

scan lines. which the. data in the interpolator buffer is retrieved and 

In order to fit this 432*512 pixel image in a space supplied over line 50 to the recorder where it is typi- 
definedby 1012.5 * 1200 spots and lines in accordance 3Q cally used to modulate the intensity of beam 42 which in 

with the present invention, we proceed as follows: In turn exposes film sheet 36. 

the CPU we calculate a ratio R corresponding to the Besides controlling the sampling rate of the data from 

ratio of the output image size in one dimension at least the buffer in the interpolator, factor f is further used to 

to the original image size in thesirne dimension. In this mo dify the speed of the film sheet movement. Factor f 

illustration, using the vertical (Y) dimension for conve- is sen t over cable 26 to a speed controller 28 connected 

mencc, this ratio is 1200 display raster lines divided by via cable 30 to motor 32. Controller 28 adjusts the speed 

512 original raster lines or, of motor 32 m a manner mverse]y p ropo J rtional t0 7 

R= 2.34375. Of course f could also be chosen as the ratio of R/I in 

which case 1/f should be substituted in the applications 

Following determination of R the two nearest inte- 40 ab 2 VCf sli11 within the of ^ mvent ion. 

gers to R are identified; These are "2" and "3". Either Returning now to the present numerical illustration, 

may be used in practicing this invention, however it is f = =I / R = I- 2 * and 1/f =0.78 125. Using this information, 

preferred that the integer which is next higher to the . ^P^S rfl te of the buffer will be increased by 1.28 

value of R be selected as experience has shown that this timcs ^ nona *} . rate at which data is supplied to the 

results in more uniform appearing output images. Here 45 recorder for writing in the X, or horizontal direction, 

this integer is 1=3. ^ . At the same time the recorder drive of the film sheet is 

The integer value is supplied through link 22 to the slowed }>y a factor of 0.78125 times the normal film 

data processing block 24 which includes an interpolator s ^ cet d " ve rate * ^ s can °e seen » 001,1 changes involve 

to which data from the original image has also been onJv smaI1 alterations in the sampling and drive rates, 

provided. Based on the value of I, the incoming data 50 which do not result in detectable artifacts or distortions, 

from memory 20 is digitally expanded using any one of * n a preferred embodiment this invention is best im- 

a number of known data expansion techniques such as plemented using software and a computer adapted to 

replication, linear interpolation, cubic spline or two operate in a system of the type shown in FIG. 4. The 

dimensional kernel convolution interpolation etc. In the following text describes the algorithms for computing 

preferred embodiment of this invention, the digital in- 55 maximally efficient format structures for the output 

terpolator employs parallel processing hardware of the image size determination, and provides the necessary 

type disclosed in co-pending Ser. application No guidance for thc preparation of appropriate software to 

07/692,655 filed concurrently herewith in the name of implement this invention in a radiological environment 

T, Byrne and assigned to the assignee of the present Should thc invention be applied in another environ- 

invention, the contents of this application being incor- 60 ment, such as publishing, including desktop publishing, 

porated herein by reference. The expanded dau is the same principles apply even though some of the 

stored in a temporary buffer which provides sufficient algorithms must be modified, particularly those that 

storage capacity for the expanded data in terms of both relate to the image distribution in the film sheet and 

horizontal and vertical expansion. resultant calculation of the optimum image output size. 

In the present numerical illustration, assuming that 65 Such modifications will be needed in order to accom- 

the digital expansion method selected is a simple repli- modate output requirements of the specific application 

cation, as each address value arrives from the stored e.g., page format, presence of more than one columns in' 

original image data to the interpolator, it is stored as a page, etc. The following process is therefore disclosed 
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here as an example of a preferred embodiment in a 
particular application; the disclosed methodology on 
the other hand may be followed to provide software 
applicable to applications other than radiology, and 
such software is considered also be within the scope of 5 
the present invention. 

In developing software to implement this invention in 
a medical radiography diagnostic environment, the 
following assumptions were made: 

1. The writing spot is substantially symmetrical in 10 
cross section; 

2. The film is scanned across its width, which is the 
shorter dimension for rectangular sheets and which 
is the X or horizontal direction referred to herein- 
above; 15 

3. The writing density (that is the number of written 
pixels per recorder spot) is, at least, equal to the 
optimum for the spot diameter. Changes in writing 
density, i.e., spots per inch, will cause greater over- 
lap of adjacent spots. The spot density is adjusted 20 
symmetrically to preserve image aspect ratio on 
magnification. If the original image pixels are not 
symmetric, the horizontal writing density is ad- 
justed to compensate; 

4. Images of different pixel composition are permitted 25 
on the same sheet, but not in the same row. In this 
case, writing density will be calculated on the basis 
of the image having the largest primitive area (pix- 
els * lines) to be displayed on a given sheet; 

5. In mixed formats for images from the same source, 30 
writing density is computed for the smaller image. 
The larger images arc displayed at the best integer 
magnification allowed by the rules defining the 
selection of borders. 

6. The transport film sheet drive will transport the 35 
film past the exposure station at a rate which re- 
mains constant during the full image display; how* 
ever, writing density may change between rows of 
images to accommodate images with asymmetric 
pixels. Otherwise, spot density will be constant for 40 
one sheet; 

7. Within the constraints of 5 and 6 t magnification by 
integer interpolation may be varied between rows 
of images. 

Following are definitions for the variables used in 45 
describing Maximally Efficient Formatting: 
FILMWID=Film width (short dimension), inches 
FILMHGT=Film height (long dimension), in 
TOPBOR = Top border, in 

BOTBOR = Bottom border, in 50 
SIDEBOR = Side border, in 
BTWNBOR= Border between images, in 
NUM WIDE = Number of images within the width of 
the film 

NUMHIGH = Number of images within the height of 55 
the film 

MXIMWID= Maximum image width, in 
MXIMHGT= Maximum image height, in 
FMTA S PR '= Format aspect ratio=MXIMWID / 

MXIMHGT 60 
IMGWID= Formatted image width, in 
IMGHGT= Formatted image height, in 
IMGPDC= Source image width, pixels 
IMGLIN= Source image height, lines 
PIXASPR= Pixel aspect ratio (usually equals one) 65 
IMGASPR^Image aspect ratio =IMGPIX • PIX- 

ASPR / IMGLIN 
SPTDENS = Spot density, spots/in 
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FIXDENS= Pixel density, pixels/in =IMGPIX / 
IMGWID 

LINDENS = Line density, lines/in = IMGLIN / 
IMGHGT 

MAGNIFY = Magnification, spots/line = SPT- 

DENS/LINDENS = (R) 
INTRPOL= MAGNIFY rounded to next highest 

integer= (I) 

DENS AD J m Density adjustment = INTRPOL / 
MAGNIFY =(0 

The first step is the determination of the desired out- 
put image size based on optimum distribution of a given 
number of source images on a preselected size film 
sheet, displayed so as to maximize use of the display 
surface area. Based on the minimum allowable border 
dimensions, the film is mapped to determine the maxi- 
mum area available for each image in a given format: 

MXIMWID~[FILMWID-l*SIDEBOR-{NUM- 
WIDE-\)*BTWNBOR)/NUMWIDE and 

MXIMHGT = [FILMHGT- TOPBOR -BOT- 
BOR-(NUMHIGH-\) 
•BTWNBOR)/NUMHIGH 

The above applies to formats having equal-size im- 
ages. For mixed formats in which there are different 
size images on the same sheet, e.g., 6/2, the calculations 
are made on the base format first, e.g., 12 images per 
sheet. The rows of base format images are arranged, 
e.g., 2 rows of 3 images, then the remaining area is filled 
with the larger images as efficiently as the sheet spot 
density will allow, by changing only the value for inte- 
ger interpolation. 

In cases where it is desired to print images for view- 
ing in horizontal (or landscape) rather than vertical (or 
portrait) orientation, the images can be rotated such that 
their pixels and lines are interchanged, and border di- 
mensions will also be different; for example, the left 
SIDEBOR dimension will equal the TOPBOR mini- 
mum in subsequent computations. 

Once the maximum dimensions of the images are 
computed, the source images selected for display are 
sized to fit in those dimensions. In typical radiological 
applications the header information for a sheet of im- 
ages contains the format selected and the dimensions of 
all the images to be displayed on that sheet Source 
image data are given as: IMGPIX, IMGLIN and PIX- 
ASPR. 

The largest source image on the sheet is used for the 
calculations: 

Largtst image=[IMGPIX • IMGLIN • 
PIXASPR)max 

The best fit is obtained when the image is magnified 
until it reaches either MXIMWID or MXIMHGT. 
Which boundary is reached first can be predicted by 
comparing the aspect ratios of the format and the image: 

If IMGASPR > FMTASPR, IMGWID** MXIMWID 
and IMGHGT** MXIMWID / IMGASPR 

If IMGASPR < FMTASPR, IMGHGT^ MXIMHGT 
and IMG WID= MXIMHGT • IMGASPR 

The magnification required to map the image pixels 
into the format area is next arrived at as. follows: 

LINDENS — IM GUN / IMGHGT 
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MAGNIFY^SPTDENS / UNDENS 

Determine INTRPOL, then compute: 

DENS A DJ = JNTR POL / MAGNIFY 5 

The DENSADJ value will be used to set the film 
transport speed and the pixel clock rate as follows: 
Under the assumption of a symmetrical spot, the hori- J0 
zontal (scanning direction) and the vertical (transport 
d irecti on) resolution are the same, i.e., the same 
SPTDENS value. If the writing density is to be in- 
creased by DENSADJ, this is accomplished by slowing 
the transport rate and increasing the pixel clock rate by 15 
the same percentage. 

Adjusted transport rate = {Nominal transport 
rate)/DENSADJ 

Adjusted pixel rate ^{nominal pixel rate) • DENSADJ 20 

The images to be printed are assembled contiguously 
in a memory large enough to hold a complete sheet. As 
described earlier in conjunction with FIG. 4, the source 
image data is concatenated with border pixels in a line 25 
buffer in a primitive state, i.e., uninterpolated, for each 
line to be written on the film. The buffer, containing 
enough lines for vertical interpolation, is integer-inter- 
polated to a second, single-line buffer which holds INT- 
RPOL * the number of primitive pixels. The second 30 
buffer is then read at the pixel clock rate in synchronism 
with the scanned spot, one line at a time. 

To calculate the number of pixels to assign to the 
borders, it is necessary to compute the primitive pixel 
density: 35 

PtXDENS-IMGPIX / JMGWJD pixels/in 
LINDENS = IMG UK / IMGHGT lines/in 

40 

From the pixel and line densities, border pixels and 
lines are simply the product of the border dimensions 
and the densities. 

For example, 

45 

TOPBOR pixels=FILMWID • PIXDENS 
TOPBOR lines= TOPBOR • LINDENS 
SIDEBOR pixels-SIDEBOR •PIXDENS. etc. 

50 

Finally there are two ways to arrange the images on 
a film sheet. Centered or grouped. In a centered ar- 
rangement the image is magnified from the center of the 
area defined by the minimum border dimensions so that 
unfilled area is symmetrically divided around the 55 
printed image. A grouped arrangement on the other 
hand uses the border between images to separate the 
images, and centers the whole group on the film. Either 
method can be accommodated as follows: 



CENTERED IMAGES: 

TOPBOR* TOPBOR + ((MXIMHGT - IMGHGTY2] 

BOTBOR# = BOTBOR + j 

BTWNBOR* = BTWNBOR + 2 • Q high 

SIDEBOR* o SIDEBOR + {(MXIMWID - IMGWID)/2} 

BTWNBOR* = BTWNBOR + 2*0 wide 

GROUPED IMAGES: 

TOPBOR* - TOPBOR + NUMHIGH • fj 



-continued 

BOTBOR* - BOTBOR + NUMHIGH • 0 

BTWNBOR* m BTWNBOR 

SIDEBOR* = SIDEBOR + NUMW1DE * {} 



In either case, the pixel and line composition of the 
adjusted borders (#) is determined by multiplying the 
computed dimensions by the appropriate spot densities. 

In most cases, the images will fill the available width 
or height on the film. However, in those cases where 
there are mixed formats or mixed images, the adjusted 
border calculations must be made in both directions to 
present the non-optimum images properly. It may be 
better aesthetically to choose centering over grouping 
for mixed images. 

Up to this point the invention has been described as 
applied to a symmetrical spot, using a square pixel. This 
is not to be construed a limitation to the claimed inven- 
tion, and asymmetrical pixels may be treated in the same 
manner with the aforementioned techniques and equa- 
tions. For example, a popular ultrasound imaging piece 
of equipment uses 768 pixels by 480 lines to define an 
image having a 4:3 aspect ratio. This means that the 
pixel aspect ratio is no longer 1.0 but 0.833. To handle 
this problem, we define R based on the assumption of 
square pixels, and modify the factor f by the pixel aspect 
ratio. R, as computed for in the Y dimension is: 

R= Display raster lines / Original Image rater lines 

fj,=sI/R. This value controls the film drive speed as 
1/f,. 

f x ={ y / Pixel aspect ratio. This value controls the 
output clock rate, that is the rate at which data is read 
out of the Interpolator. 

Thus, in cases where the pixel aspect ratio is different 
than 1.0, calculations as shown in the above example, 
must be made in both the X and Y dimensions. 

Those having the benefit of the teachings of our in- 
vention as set forth above may attempt various modifi- 
cations and applications not expressly set forth herein; 
thus* for example, the same principles and methods 
taught hereinabove may be equally advantageously 
applied to a multi-plane representation as is encountered 
in 3 or 4 color systems, and still remain within the scope 
of our invention. The invention is not limited to the 
specific embodiments described above, but rather is 
defined by the appended claims and their fair equiva- 
lents. 

We claim: 

1. A process for fitting a source image comprised of 
raster image data representing an arbitrary number of 
pixels and an arbitrary aspect ratio, in a specified raster 
output image having a desired size and aspect ratio 
comprising: 

(a) determining a desired output image size in at least 
one dimension, 

(b) determining a source image size in the same at 
least one dimension, and comparing said source 
image size to said output image size to obtain a ratio 
R in that at least one dimension; 

(c) calculating a first integer I adjacent to the ratio R; 

(d) digitally expanding the source image by expand- 
ing the source image data by a factor equal to the 
calculated integer I to generate an expanded output 
image data; 

(e) calculating the quotient of the calculated integer I 
to the ratio R to obtain a sampling factor f where 
f=I/R; 
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(0 sampling the expanded output data at a sampling 
rate controlled by the sampling factor f, and 

(g) using the sampled output data to drive an output 
display device. 

2. The process in accordance with . claim 1 wherein 5 
the the integer I is the next higher integer to the ratio R. 

3. The process in accordance with claim 1 or 2 
wherein the raster scan system scans a surface at a scan 
rate and wherein the factor f=I/R is also used to con- 
trol the raster scan rate. 10 

4. The process in accordance with claims 1 or 2 
wherein the output data is used to modulate a writing 
beam intensity and wherein the intensity of the beam is 
also adjusted to compensate for a change in exposure 
needed to compensate for a corresponding change in 15 
the line spacing. 

5. The process in accordance with claim 1, or 2, 
wherein the digital expanding step is performed using 
any one of the following methods: 

(a) Replication, or 20 

(b) Bilinear interpolation, or 

(c) Cubic-spline interpolation, or 

(c) Two dimensional kernel convolution interpola- 
tion. 

6. The process of claim 1 or 2, wherein the desired 25 
output image size is determined based on optimum dis- 
tribution of a given number of source images on a prese- 
lected size film sheet. 

7. The process of claim 1 or 2, wherein the source 
images represent medical diagnostic images. 30 

8. The process of claim 1, wherein the source image 
size in step (b) is determined in each dimension used, 
and the source image size is compared to the output 
image size sought, to obtain image ratios in each dimen- 
sion used, and wherein the integer I calculated in step 35 
(c) is the first higher integer adjacent to the larger of the 
two image ratios determined, said larger ratio being 
identified as R. 

9. A process in accordance with claim 1 wherein the 
source image pixels have an aspect ratio other than one, 40 
and the source and output images comprise an arbitrary 
number of raster lines each, and wherein, 

R= ratio of the display raster lines to the original 
raster lines. 

45 
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10. A process in accordance with claim 9, wherein 
the factor f is used to control the raster scan rate of the 
display device, and the ratio of the factor f to the pixel 
aspect ratio is used to control the sampling rate of the 
expanded output data. 

11. Apparatus for use with a raster scan type image 
reproduction recorder having control means to adjust a 
scanning rate in at least one dimension, the equipment 
adapted for accessing a source image comprised of: 
source image data represented by an arbitrary number 
of pixels and an arbitrary aspect ratio, the data residing 
in a storage means, the equipment further adapted for 
generating data for reproducing in said recorder said 
image in a specified output image format having a de- 
sired size and aspect ratio comprising: 

(a) a central data processing means including means 
for accessing a preselected image data in a storage 
medium; 

(b) means also in said central processing means for 
calculating: 

i. a ratio R where R= output image size/source 
image size; 

ii. a first integer I adjacent, next lower or next 
higher, to the ratio R; and 

iii. a quotient of the ratio R to the calculated inte- 
ger I to obtain a sampling factor f where f=I/R; 

(c) means controlled by said central processing 
means, for retrieving said preselected image data 
and for delivering said data to a data interpolating 
and storing means adapted to digitally expand the 
source image data by a factor equal to the calcu- 
lated integer I to generate expanded output image 
data and to store at least a portion of such data; 

(d) means associated with the central data processing 
means to apply the sampling factor f to a sampling 
rate control means to adjust the sampling rate of 
the expanded output data generated by the interpo- 
lator, and to supply said data to the recorder; and 
to also apply the sampling factor f to the control 
means for adjusting the vertical scanning rate of 
the recorder. 

12. The apparatus of claim 11, wherein the first inte- 
ger I is the next higher integer to the ratio R. 
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